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We  study  the  singular  Cauchy  problem 


( x,  t )  e  R  X  [0, T] ,  T  >  0 


1 Ad 


,  x  e  r  . 

max  (0,  Or  the  initial  datum  g  is  smooth 


v  =  $(v)  , 

t  xx 

v(x,0)  -  g(x) 

The  constitutive  function  $(0 

(bounded)  on  R\{0}  and  satisfies 

g(0)  »  0j  xg(x)  >  0,  x  e  Rr  g*(0+)  •  g' (0~)  f  0  , 

where  the  superscript  "+"  denotes  the  limit  from  the  right  (left). 

show  that  the  free  boundary  s  given  by  v(s(t)+,t)  =*  0  satisfies 


We 


s(t)  -  -ic/t  +  o(/t)  (t  +  0+)  , 

where  <  >  0  is  a  monotone  function  of  p  *  g' (0+)/g’ (0  ),  implicitly 
defined  by  the  equation 


(*) 


P 


k  >  0 


This  generalises  an  earlier  analysis  [7]  for  the  special  case  of  smooth  data 
g  in  which  g' (0+)  -  g'(0“)  +  0,  and  p  -  1.  In  this  case  the  numerical 
value  <  -  .903446...  as  computed  from  <*)  is  consistent  with  our  previous 
result . 

AMS  (MOS)  Subject  Classifications:  35K055,  35K65,  45G05 

Key  Words:  Cauchy  problem,  parabolic,  nonlinear,  regularity  of  free  boundary, 
self  similar  solutions,  nonlinear  singular  integral  equation. 
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SIGNIFICANCE  AND  EXPLANATION 

^  The  Cauchy  problem  stated  in  the  Abstract  is  similar  to  the  well-known 
one  phase  Stefan  problem  (in  one  space  dimension).  In  the  latter  one  would 
assume  g(x)  =  -1  for  x  <  0,  as  well  as  g(x)  >  0  for  x  >  0,  so  that 
g  would  have  a  jump  discontinuity  at  x  =  0.  Our  assumptions  on  the  initial 
data  g  yield  a  different  behavior  of  the  solution  v  and  of  the  resulting 
free  boundary.  Indeed,  the  free  boundary  is  not  (infinitely)  differentiable 
at  t  ■  0,  contrary  to  the  situation  for  the  classical  Stefan  problem. 

This  problem  also  serves  as  a  prototype  of  nonlinear  parabolic  problems 


r 

which  arise  as  monotone  ^convexifications  of  nonlinear  diffusion  equations 


o' 


with  nonmonotone  constitutive  functions  $  (see  [5],  [6]).  That  analysis 
shows  the  existence  of  infinitely  many  solutions  v  of  the  nonmonotone 
problem  each  having  v  bounded,  but  oscillating  more  and  more  rapidly  as 


oCt) 


A 

Thus  each  solution  v  exhibits  phase  changes.  Numerical  experiments 


t  ♦  sf . 

further  suggest  the  conjecture  that  the  "physically  correct^  solution  of  the 
nonmonotone  problem  is  the  one  which  for  t  >  0  sufficiently  large  approaches 
the  unique  solution  of  the  appropriately  related  convexified  monotone 
problem.  This  paper  is  another  step  towards  the  understanding  of  this 

intriguing  phenomenon.  Our  earlier  analysis  of  the  Cauchy  problem  stated  in 

A' 

the  abstract  was  for  smooth  data  g  with  g’ (0+)  ■  g’ (0  )  >  0,  [7).  The 

present  analysis  shows  how  the  free  boundary  s(t)  depends  near  t  *  0  on 
more  general  data  for  which  p  ■  g’ (0+)/g’ (0  )  need  no  longer  be  1.  We 
present  two  approaches:  (i)  a  generalization  of  the  integral  equation  for  the 
free  boundary  studied  in  [7],  (ii)\a  preliminary  analysis  via  self  similar 


solutions  which  will  be  exploited  in  .  [8]. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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10]  where  g(x)  “  -1  for  x  <  0.  However,  the  assumption  (2)  yields  a 
different  behavior  of  the  free  boundary  and  of  the  solution  v.  In  fact,  for 
smooth  initial  data  (g* (0+)  “  g' (0  )),  vx  is  not  continuous  at  (x,t)  “ 
(0,0),  whereas  solutions  to  the  Stefan  problem  are  smooth  on  the  set 
{ ( x , t )  :  x  >  s(t),  t  e  (0, T] } . 
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For  a  sufficiently  regular  solution,  problem  (1)  -  (2)  is  equivalent  to 


the  free  boundary  problem 

^  v  =  v  ,  s(t)  <  x  <  «• 

t  XX 

\  v(s(t) ,t)  =  0 
I  v(x,0 )  =  g(x) 


f  g(s(t))s'(t)  =  v  (s(t),t),  t  >  0  , 

s(0 )  -  0  . 


One  easily  verifies  that  a  solution  of  (2)  -  (4)  which  is  extended  by 
(5)  v(x,t)  -  g(x)  ,  -»  <  x  <  s(t)  , 

is  a  weak  solution  of  problem  (1),  (2). 

In  [7]  we  studied  the  equivalent  integrated  version  of  problem  (1)  in  the 


form 

(!') 


u  «  $(u  ) 
t  XX 


u(x,0  5  *  f(x)  :=  /*  g(y)dy 


Assuming  sufficient  regularity,  the  solutions  of  problems  (1)  and  (1‘)  are 
related  by 

u(x,t)  =  J*  v(y,t)dy  +  jjj  $(v)x(0,T)dT 
The  boundary  condition  (4)  can  be  also  easily  derived  from  (1'): 
continuity  of  u  across  the  free  boundary  s  implies 

f  ( s(  t) )  =  u(  s(  t) ,  t) 

Differentiating  with  respect  to  t  and  using  ux(s(t)+,t)  •  v(s(t)+,t)  •  0 
and  f’(x)  *  g(x)  we  obtain  (4). 


We  assumed  in  [7]  that  the  initial  data  f  for  the  Cauchy  problem  (1*) 


For  smooth  initial  data  (i.e.  g'fO-)  -  g' (0+))  we  have  p  ■  1  and  x  m 
0.9034...  which  is  consistent  with  our  result  in  [7]. 

A  complete  proof  of  the  ‘Rieorem  will  be  included  in  a  forthcoming  joint 
paper  with  J.  Vazquez  [8].  It  is  based  on  comparison  arguments  using  self 
similar  solutions  of  problem  ( 1 ) .  More  generally,  we  shall  consider  initial 
data  of  the  form 


g(x)  -  < 


x  g+(x) 


m 


x  >  0  , 


^ -|x|  g_(x)  ,  x  <  o  , 

with  9^(0*)  >  0  and  determine  the  local  behavior  of  s  in  all  cases. 

In  this  note  we  prove  the  Theorem  by  extending  the  method  used  in  [7]. 

This  approach,  compared  to  that  using  similar  solutions,  yields  additional 

regularity  for  the  free  boundary,  namely  for  any  o  e  (0,^)  we  have 

(8)  j-  [g(s(t))s*(t)]  -  Oft®"1),  as  t  0+  . 

at 

However,  for  technical  reasons,  we  have  to  assume  that  p  <  1. 

We  begin  by  constructing  in  Section  2  self  similar  solutions  to  problem 
(2)  -  (4)  and  derive  the  relation  (7).  The  behavior  of  these  solutions  is 
typical  for  the  general  case  treated  in  Section  3.  We  prove  (6)  and  (8)  by 
solving  a  nonlinear  integral  equation  for  the  function  r(t)  g(s( t) )s' (t) . 


2.  Self  similar  solutions . 

For  the  special  initial  data 


(9) 


it  will  be  shown  below  that  problem  (2)  - 


x  >  0  , 

x  <  0  , 

(4)  has  self  similar  solutions  of 


the  form 


(10) 


I 


v<x,t)  -  t^*(xt”  1/2  )  . 

Substituting  (10)  into  equations  (3)  and  (4)  and  assuming  that  s(t)  * 
-x/t  we  see  that  must  satisfy  the  ordinary  differential  equation 

(11)  ♦"( C)  +  \  C  ♦  '(C)  -  J  WO  -  0,  C  >  -x  , 
with  the  boundary  conditions 

(12)  ♦<— *c)  -  0,  ♦'(-x)  -  p_  j  k2 

(13)  lim  c“Vo  -  p.  • 

C— 

Note,  that  (13)  is  equivalent  to  the  initial  condition  v(x,0)  m  p+x. 
Observing  that  c^C  is  a  particular  solution  of  (11)  for  any  constant  c^, 
one  easily  finds  that  the  general  solution  of  (11)  is 

♦(C)  -  c^C  +  c2(2e“*  /A  +  C  J*  e"*  /4dy)  • 

From  (12)  we  obtain 


.X  ,  X  X  /4  r* 

?.(r +  r  *  K 


-y  /4 


dy) 


3  2  ,A 

xx/ 4 

p.  T  * 


Therefore,  by  (13),  p+,  p_  and  x  must  satisfy 

C1  +  C2^W  “  p+  * 

which  is  the  relation  (7).  This  establishes  the  existence  of  self  similar 
solutions  for  the  model  initial  data  (9).  Since  s(t)  -  -x/t  the  assertions 
(6)  and  (8)  of  the  Theorem  are  trivially  valid  in  this  case. 


? 


i 


I 


L 


3.  Proof  of  the  Theorem. 


Existence  and  uniqueness  of  weak  solution  for  problem  (1)  follow  from 

nonlinear  semigroup  theory  [2].  It  la  therefore  sufficient  to  show  existence 

of  a  solution  (v,s)  for  problem  (2)  -  (4)  and  prove  the  assertions  (6)  -  (8) 

for  the  free  boundary  s.  Making  the  restrictive  assumption  p  <  1  this  can 

be  done  by  extending  the  method  in  [7], 

1  x2 

Let  r (x,t )  :*  -  exp(-  — )  denote  the  fundamental  solution  of  the 

✓  m  4t 

heat  equation.  Assume  that  (v,s)  is  a  solution  of  problem  (2)  -  (4)  for 

which  the  function 

(14)  r(T)  :=  v  <s(t),t) 

x 

is  bounded  and  continuous  on  ( 0, t] .  Integrating  Green's  identity 

(r<x-£,t-T)w5(?,T)  -  rx(x-£,t-T)w(C,T))c  -  (r(x-£,t-T)w(£,T))T  =  o 

over  the  domain  {(£,t)  :  s(t)  <  5  <  •,  T  6  (0,t)},  following  [7]  one  can 
derive  an  integral  equation  for  r: 

(15)  r(t)  =2  /"  r(s(t)-5,t)g' (C)d?  -2  /*  I*  (  s( t )-s( t) ,t-r)r ( x)dT  . 

With  the  abbreviation 


A(s,t,T) 


s(t)-s(tT) 

- W 

2(t-tT)/2 


this  equation  can  be  rewritten  as 

r(t)  «  —  /“  exp(-  7  (S^  -  5>2)g’U  /t)dC 
✓x  0  4  /t 

(16)  +  —  fl  T -  exp(-A(s,t,T)2)r( tT)dr 

0  1_T 


(Fr ) (t)  +  (Kr)(t)  . 


The  Theorem  is  a  consequence  of  the  following  existence  result  for  (16). 
Proposition.  Assume  that 

p  -  g'(0+)/g'(0* )  — .  p+/p_  <  1  . 


-6— 


Let  s  and  r  _in  (16)  be  related  by 

(17)  r(t)  =  g(s(t) )s* (t)  . 

Then ,  for  any  a  e  (0,1/2),  there  exists  T  >  0  such  that  the  integral 


equation  (16)  has  a  solution  T  for  which 

(18)  sup  t1  ° | f *  < t ) I  <  •  . 

te(0,T) 

We  postpone  the  proof  of  this  proposition  and  complete  the  proof  of  the 
Theorem. 


Let  r  be  a  solution  of  (16)  for  which  (18)  holds.  Then,  it  is  easy  to 

see  that  s,  given  by  (17),  satisfies 

—  1/.  +o 

(19)  s(t)  +  K/t  *  0(t  2  )  , 

1/, 

with  k  =  (2r(0)/p  )  2.  To  verify  (7)  we  pass  to  the  limit  t\  0  in  (16). 


Since 


we  obtain 


A(s,0,T) 


K  1-/t 


9 


(20) 


“  /"  exp(-  ^(«+C)2)p+  d? 

/if 


1  rl  1C  1  .1C 

-  fo  2  exp(-  r 


1  — /  T 


/1-t<hVt 


1+/t 


±>P_  o  dx 


Substituting 

get 


in  the  second  integral  and  solving  for 


P 


P+/P-  we 


(20' ) 


(—  +  —  —  f1  lat 

'2  4  Jo  : 

1+y 


Mr. 

'  u  K 


^-y2/4 


dy 


1 


* 


-7- 


A  simple  calculation  shows  that  this  is  equivalent  to  (7):  Subtracting  the 


right  hand  sides  of  equations  (7)  and  (20')  and  multiplying  by  /*  x 
gives 

,1  .  1  x2/4  r<*  -y2/4  /•  -y2/4  /*  1  1  -y2  -x2y2/4 

(2*  4  e  /»k  6  dy)  *  Ac  e  dy  ‘  ^  +  ?  / 0  2  6  dy) 

i-*y 

Dsing  e^  ^4dy  =  /*  -  /*  e^  ^4dy,  the  last  expression  simplifies  to 

1  ie2/4r  rl  4 _  -x2 (y2+1  >/4^  .  ,  .  fx  -y2/4^  ,2,. 

4  e  t-  /0  - j  e  dy  +  (x  -  <J0  e^  dy)  )]  . 

14y 

The  term  in  square  brackets  vanishes  for  x  ■  0.  Since  also 


—  [•••] 
dx 


2  /;  ,-y2/4< 


this  finishes  the  argument. 

It  remains  to  prove  the  Proposition.  To  this  end  we  define  for 

oe  (0,1/2)  the  seminorm 

I r |  :*  sup  t1  a|r'(t)| 

1  te(0,T) 

Further,  we  denote  by  c  a  generic  constant  which  may  depend  on  g,  o,  x  and 
monotonely  on  T  and  | r 1 T . 

2 

If  r{0)  ■  PJ2~  with  K  9i^en  by  (7)  and  s  is  related  to  r  by  (17), 
we  claim  that  the  following  estimates  hold: 

(21)  |FrlT  <  c(T)  +  (Cl(x,a)  +  c(T))  |rlT 

(22)  |Xr|T  4  (c2(x,o)  +  c(T ) )  | r 1 T 


■8- 


with 


c^KfO)  -  “  E-^-  exp{-<2/4) 


c2(ic,a)  - 

— “  <  fl - ~ - —  exp(-  j-  1-^)dT  + 

2/w  /l-T  (1+/t  1+A 


(2+2a)(2+4a) 


fl  1-t1/2+* 
'°  (H)3/2 


1  «  K2  ,  »«2  1-/t.  , 

M  “  T - “I  exp(-  - - -)dT  . 

1+/t  4  1+/t 


Here,  T  c(T)  denotes  a  continuous  function  with  c(0)  -  0.  Note,  that  the 
constants  c^  and  c2  do  not  depend  on  T  and  hence  are  not  small  for 
small  T.  this  reflects  the  fact  that  the  operators  P  and  K  are  not 
compact.  Our  existence  proof  relies  on  the  fact  that 
(23)  w  :*  c^x,'-)  +  c2(x,^)  <  1 

which,  however,  is  valid  only  for  it  <  1.05...,  or  equivalently  p  <  1.59...  . 
this  explains  the  artificial  restriction  on  p  in  the  Proposition. 

Assuming  (23),  we  can  proceed  in  a  standard  manner.  We  iterate  the 
integral  equation  (16)  in  the  form 


V*3  5  P_ 

(24) 

rn+1  "  ^n  +  "n  '  n  e  W  ' 

with  p  related  to  x  by  (7).  By  (20)  r  (0)  -  r,  for  all  u.  Moreover, 

n  i 

if  we  choose  a  close  to  1/2  and  T  >  0  so  small  that 

(c,|(ic,a)  +  c(T ) )  +  (c2(x,a)  +  c(T))  -s  w’  <  1  , 

the  estimates  (21)  and  (22)  imply 


It  follows  that 


I'n+i1!  <  C<T)  +  W|rnlT 


I rn  lT  *  c(T)  /  H-**’  )>  «  e 


and  we  can  select  a  subsequence  which  converges  to  a  solution  r  of  (16) 


in  C[0,T)  which  satisfies  (18). 


9 


Bie  proof  of  the  estimates  (21),  (22)  is  fairly  technical  and  we  refer  to 


[7],  since  only  minor  modifications  are  necessary  to  incorporate  the  case 
g* (0+)  y  g*  (0~).  It  is  not  quite  satisfactory  that  one  has  to  compute  the 
constants  c1  and  c2  explicitly.  It  may  perhaps  be  possible  to  avoid  this 
by  estimating  the  nonlinear  operators  F  ar.d  K  in  a  suitable  weighted 
norm.  In  fact  we  conjecture  that  a  stronger  result  can  be  proved  namely  that 
for  smooth  initial  data  g,  r  and  s  are  smooth  function  of  /t  at  t  *  0. 


-10- 
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t  XX 

v(x,0)  ■  g(x)  ,  x  6  R  . 

The  constitutive  function  <M€)  -  max(0,t)j  the  initial  datum  g  is  smooth 
(bounded)  on  R\{0>  and  satisfies 

i0)  -  0*  xq (x)  >.0,  x  6  R>  g'  (0*)  •  g1  (0  )  0  0  , _ 
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ABSTRACT  (continued) 


where  the  superscipt  "+"  ("-••)  denotes  the  limit  from  the  right  (left).  We 
show  that  the  free  boundary  s  given  by  v(s(t)  +  ,t)  *  0  satisfies 

s(t)  ■  -<</t  +  °(Jt)  (t  -*•  o+)  , 

where  tc  >  0  is  a  monotone  function  of  p  =  g'  (0+)/g'  (0~),  implicitly 


defined  by  the  equation 


2  3  2  ..  .2  .. 

K  tc  <  /4  /«  -t  /4 
-r-  +  ~T  e  I  e  dt,K>0 
2  4  '  -< 


This  generalizes  an  earlier  analysis  [  ]  for  the  special  case  of  smooth  data 

g  in  which  g* (0+)  =  g' (0  )  ?  0,  and  p  »  1.  In  this  case  the  numerical 
value  tc  =  .903446...  as  computed  from  (*)  is  consistent  with  our  previous 
result. 
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